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Abstract. In this paper, indirect collocation approach based on com¬ 
pactly supported radial basis function (CSRBF) is applied for solving 
Volterra’s population model. The method reduces the solution of this 
problem to the solution of a system of algebraic equations. Volterra’s 
model is a non-linear integro-differential equation where the integral 
term represents the effect of toxin. To solve the problem, we use the 
well-known CSRBF: Wendlandz^^. Numerical results and residual norm 
(||R(t)|P) show good accuracy and rate of convergence. 
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1. Introduction 


The volterra’s model for population growth of a species within a closed 
system is given in [1, 2] as 


( 1 . 1 ) 


dp 

— = ap-bp^ -cp J p{x)dx, p{0) = po, 


where a > 0 is the birth rate coefficient, 6 > 0 is the crowding coefficient 
and c > 0 is the toxicity coefficient, the coefficient c indicates the essential 
behaviour of the population evolution before its level falls to zero in the long 
term. pQ is the initial population and p = p{t) denotes the population at 
time t. This model includes the well-known terms of a logistic equation, and 
in addition it, includes an integral term cp fQp(x)dx that characterizes the 
accumulated toxicity produced since time zero [2, 3]. 

We apply scale time and population by introducing the non-dimensional 
variables 


( 1 . 2 ) 


tc pb 

t = u = 
b a 


to obtain the non-dimensional problem 


(1.3) 


du 

dt 


= u — u — u 


u{x)dx, u{0) = uo, 
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where u{t) is the scaled population of identical individuals at time t and 
K = ^ is a prescribed non-dimensional parameter. The only equilibrium 
solution of Eq. (1-3) is the trivial solution u{t) = 0 and the analytical 
solution [2] 

1 f*- r 

(1.4) u(t) = uo exp(— / [1 — M(r) — / u{x)dx]dT) 

^ Jo Jo 

shows that u{t) > 0 for all t if uq > 0. 

The solution of Eq. (1.1) has been of considerable concern. Although 
a closed form solution has been achieved in [1, 5], it was formally shown 
that the closed form solution cannot lead to any insight into the behaviour 
of the population evolution [1]. Some approximate and numerical solutions 
for Volterra’s Population Model have been reported, the successive approx¬ 
imations method was suggested for the solution of Eq. (1.3), but was not 
implemented. In this case, the solution u{t) has a smaller amplitude com¬ 
pared to the amplitude of u{t) for the case k <C 1. 

In [2], several numerical algorithms namely the Euler method, the mod¬ 
ified Euler method, the classical fourth-order Runge-Kutta method and 
Runge-Kutta-Fehlberg method for the solution of Eq. (1.3) are obtained. 
Moreover, a phase-plane analysis is implemented. In [2], the numerical re¬ 
sults are correlated to give insight on the problem and its solution without 
using perturbation techniques. However, the performance of of the tra¬ 
ditional numerical techniques is well-known in that it using provides grid 
points only, and in addition, it requires a large amounts of calculations. 

In [3], the series solution method and the decomposition method are im¬ 
plemented independently to Eq. (1.3) and to a related non-linear ordinary 
differential equation. Furthermore, the Padeapproximations are used in the 
analysis to capture the essential behaviour of the populations u{t) of iden¬ 
tical individuals and approximation of Umax and exact value of Umax for 
different n were compared. Small [5] solved the Volterras Population Model 
by the singular perturbation method. This author scaled out the parameters 
of Eq. (1.1) as much as possible and considered four different ways to do 
this. He considered two cases k = ^ small and u = -^ large. 

It is shown in [5] that for the case k <C 1, where populations are weakly 
sensitive to toxins, a rapid rise occurs along the logistic curve that will 
reach a peak and then is followed by a slow exponential decay. And, for 
larger, where populations are strongly sensitive to toxins, the solutions are 
proportional to sech?{t). 

In [6] Adomian decomposition method and Sinc-Galerkin method were 
compared for the solution of the same integral equation in this paper. This 
showed that Adomian decomposition method is more efficient and easier to 
use for the solution of Volterras Population Model. 

In [9] the approach is based upon composite spectral functions approx¬ 
imations. The properties of composite spectral functions consisting of few 
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terms of orthogonal functions utilized to reduce the solution of the Volterras 
model to the solution of a system of algebraic equations. 

In [7] Rational Chebyshev and Hermite functions collocation approach 
were compared for the solution of Volterras Population Model growth model 
of a species within a closed system. They reduced the solution of this prob¬ 
lem to the solution of a system of algebraic equations. 

In [15] applied two common collocation approaches based on radial basis 
functions to solve Volterras Population Model. 

In [18] a numerical method based on Hybrid function approximations 
was proposed to solve Volterras Population Model. These Hybrid functions 
consist of block-pulse and Lagrange-interpolating polynomials. 

Also, in [19] Volterras population growth model of a species within a closed 
system is approximated by collocation method based on two orthogonal 
functions, Sine and Rational Legendre functions. Momani et al. [20] and 
Xu [21] used a numerical and analytical algorithm for approximate solutions 
of a fractional population growth model, respectively. The first algorithm is 
based on Adomian decomposition method (ADM) with Pad approximants 
and the second algorithm is based on homotopy analysis method (HAM). 

2. ICSRBF METHOD 

2.1. CSRBF. Many problems in science and engineering arise in infinite 
and semi-infinite domains. Different numerical methods have been proposed 
for solving problems on various domains such as FEM[13, 14], FDM[4, 13] 
and Spectral[17, 7] methods and meshfree method[26, 16]. The use of the 
RBF is an one of the popular meshfree method for solving the differential 
equations [10, 12]. For many years the global radial basis functions such 
as Gaussian, Multi quadric. Thin plate spline, Inverse multiqudric and etc 
was used [11, 8, 27]. These functions are globally supported and generate a 
system of equations with ill-condition full matrix. To convert the ill-condition 
matrix to a well-condition matrix, CSRBFs can be used instead of global 
RBFs. CSRBFs can convert the global scheme into a local one with banded 
matrices. Which makes the RBF method more feasible for solving large-scale 
problem [22]. 

Wendland’s functions. The most popular family of CSRBF are Wend- 
land functions. This function introduced by Holger Wendland in 1995 [23]. 
Wendland starts with the truncated power function <^i(r) = (1 — r)[)_ which 
be strictly positive definite and radial on for / > [|J -|- 1 , and then he 
walks through dimension by repeatedly applying the operator I. 

Definition [24] with <^i(r) = (1 — r)[)_ we define 

(2-1) 

it turns out that the functions (j)s,k are all supported on [0,1]. 

Theorem 1 [24] The function are strictly positive definite (SPD) and 
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radial on and are of the form 


(t>s,k{r) 


Ps,k{r) rE [0,1], 

0 r > 1, 


with a univariate polynomial ps,k of degree [|J + 3A: + 1. Moreover, cfis^k £ 
C‘^^{R) are unique up to a constant factor, and the polynomial degree is 
minimal for given space dimension s and smoothness 2k [24]. Wendland 
gave recursive formulas for the functions (pg^k for all s, k. We instead list the 
explicit formulas of [25]. 

Theorem 2 [24] The function (j)s^k: k = Q, 1, 2, 3, have form 

(/)s,o = (1 - 

05,1 = (1 ~ ry+^[{l + l)r + 1], 

^s ,2 = (1 - ry+'^iil^ + 4/ + 3)r2 + {31 + 6)r + 3], 

0, 3 = (1 - r)[+3[(/3 + 9^2 + 231 + 15)r3 + (6/^ + 361 + 

+(15/ + 45)r + 15], 


where I = [|J + k + 1, and the symbol = denotes equality up to a multi¬ 
plicative positive constant. 

The case k = 0 follows directly from the definition, application of the defi¬ 
nition for the case k = 1 yields 


4>s,iir) = = / t(f)i{t)dt 

Jr 

poo p\ 

= i(l — = / t{l — t^dt 

Jr Jr 

il-ry+^[il + l)r + l], 


{I -|- 1)(/ -|- 2) 


where the compact support of (j)i reduces the improper integral to a definite 
integral which can be evaluated using integration by parts. The other two 
cases are obtained similarly by repeated application of I. [24] We showed the 
most of wendland functions in table 1 . 


Table 1. Wendland’s compactly supported radial function for various 
choices of k and s=3. 


4^s,k 

smoothness 

SPD 

<t>3,oir) = (1 - r)'^ 



4>3,l{r) = (1 - r)^(4r + 1) 

C'^ 

R^ 

4>3,2{r) = (1 - r)!i(35r^ + 18r + 3) 


W~ 

03 , 3 (r) = (1 — r)^(32r^ + 25r^ + 8 r + 1 ) 



</)3 4 (r) = (1 - r)l^“(429r"‘ + 450r^ + 210r^ + 50r + 5) 


R^ 

<^3 5 (r) = (1 - r)b^(2048r5 + 2697r'‘ + 1644r^ + 566r^ + 108r + 9) 


W~ 
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2.2. Interpolation by CSRBFs. The one-dimensional function y{x) to 
be interpolated or approximated can be represented by an CSRBF as 

N 

y{x) ^ yn{x) = 

i=l 


where 


Mx) = 



= [(j)l{x), Mx), ■ ■ ■ , (i>Nix)], 

s = [6)6) • • • 


( 2 . 2 ) 

X is the input, is the local support domain and 6s are the set of coefficients 
to be determined. By using the local support domain, we mapped the 
domain of problem to CSRBF local domain. By choosing N interpolate 
nodes {xj, j = 1,2, - ■ ■ ,N) in domain: 

N 

yj = '^Ci4>i{xj)ij = 1)2, ••• ,A^). 

i=l 

To summarize the discussion on the coefficients matrix, we define 
(2.3) = y, 

where : 

y = [yi)2/2)--- ,yNf, 

A = [$'^(xi), ^>'^(X2), • • • , ^^{xN)f 
/ (t>l{xi) 4>2ixi) ••• (pNixi)\ 

4>i{x 2) 4>2ix2) ••• (pNix2) 

\4>i{xn) 4>2ixN) ••• 4>n{xn)/ 

Note that 4>i{xj) = <?i>(-^^—^-i^), by solving the system (2.3), the unknown 

coefficients 6 will be achieved. 


2.3. ICSRBF method. In the indirect method, the formulation of the 
problem starts with the decomposition of the highest order derivative under 
consideration into CRBF. The obtained derivative expression is then inte¬ 
grated to yield expressions for lower order derivatives and finally for the 
original function itself. 

We approximate ^ for solving the model by ICSRBF: 


N 

UN{t) = ^66(6 = ^^(i)S, 
2=1 


(2.4) 


du 

dt 








6 


PARAND AND HEMAMI 



Figure 1. Plot of approximate solutions of Eq. (1.3) for 
uq = 0.1 and n = 0.02, 0.04, 0.1, 0.2,0.5. 


by using integral operator f{x)dx we have 

/ du ~ / U]sf{v)dv = 

Jo Jo 

(2.5) u{t) = 44>^(t) “ + uo, 

( 2 . 6 ) I^u = +uot. 

Now, to obtain we define the residual functions by substituting Eqs. 

(1.4)-(2.6) in Eq. (1.3) 


(2.7) R{t) = + lJ{t)E-uo- 

The set of equations for obtaining the coefficients come from equal¬ 
izing Eq. (2.6) to zero at N interpolate nodes fro™ tj = j = 

1, 2, • • • ,N where L is a last interpolate node and p is a arbitrary parameter. 

(2.8) R{tj) = 0, j = 1,2,-■■ ,N. 


3. Application 

We applied the method presented in this paper to examine the mathe¬ 
matical structure of u{t). Table (2) shows the maximum of u{x) for some k 
and uq = 0.1 by using in comparison with exact solution and ADM solution 
by Wazwaz [3]. The resulting graph of Eq. (1.3) is shown in Eig. (1). 

The local support domain and arbitrary parameter p must be specihed 
by the user.An important unsolved problem is to find a approach to deter¬ 
mine the optimal size of r^[22]. The accuracy of these CSRBE depends on 
the choice of and p. By the meaning of residual function in case of Eq. 
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Table 2. A comparison of ADM[3] and the present method with exact 
values for Umax- 


K 

'^max 


P 

N 

ICSRBF 

ADM 

0.02 

0.9234271 

1 

1.766000 

15 

0.92342716 

0.9234270 

0.04 

0.8737199 

1 

1.780000 

18 

0.8737193 

0.8612401 

0.1 

0.7697414 

1 

1.811000 

18 

0.7697414 

0.7651130 

0.2 

0.6590503 

2 

1.032770 

18 

0.6590493 

0.6579123 

0.5 

0.4851902 

2 

1.114035 

27 

0.4851903 

0.4852823 


Table 3. Minimum value of ||-Res|p which is obtained with 
Tuj and p for ICSRBF. 


K 

llResf-ICSRBF 

0.5 

2.11e-08 

0.2 

2.87e-07 

0.1 

1.34e-07 

0.04 

7.86e - 05 

0.02 

4.41e - 05 


(2.6), we try to minimize ||R(t)|p by choosing good and p parameters. 
We define ||i?(t)|p as 



were 


COi = 


L 


{I - S^^{£Pm+lis)\s=s,y- 


j = ,m,Pm+i{sj) = 0, j = 0,1, ■ 


Pm,+i{x) is (m + l)th-order Legendre polynomial. Table (3) show the mini¬ 
mum of ||i?(t)|p which is obtained with local support domain and arbi¬ 
trary parameter p . 


4. Conclusion 

A method has been presented for solving Volterras Population Model 
which is an integro-ordinary differential equation, based on the compactly 
supported radial basis functions approximation. In this work, we applied 
two common ICSRBF methods on the Volterras Population Model without 
converting it to an ordinary differential equation. We used Wendland 3^5 
function. This function are proposed to provide an effective but simple 
way to improve the convergence of the solution by collocation method. As 
appeared from the Figures, we have shown that, when the constant = -^ 
is small, this type of population is relatively insensitive to toxins, and when 
c = ab is large, population of this type are extremely sensitive to toxins. 
Additionally, through the comparison with ADM, we have showed that the 
ICSRBF approach have good reliability and efficiency. 
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